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Abstract 

Using high-resolution direct numerical simulations, the height and Reynolds num- 
ber dependence of higher-order statistics of the energy dissipation rate and local 
enstrophy are examined in incompressible, fully-developed turbulent channel flow. 
The statistics are studied over a range of wall distances, spanning the viscous sub- 
layer to the channel flow centerline, for friction Reynolds numbers Re T = 180 and 
-Re T = 381. The high resolution of the simulations allows dissipation and enstrophy 
moments up to fourth order to be calculated. These moments show a dependence 
on wall distance, and Reynolds number effects are observed at the edge of the log- 
arithmic layer. Conditional analyses based on locations of intense rotation are also 
carried out in order to determine the contribution of vortical structures to the dis- 
sipation and enstrophy moments. Our analysis shows that, for the simulation at the 
larger Reynolds number, small-scale fluctuations of both dissipation and enstrophy 
become relatively constant for z + > 100. 

Key words: Turbulent channel flows, turbulent shear flows, energy dissipation 
rate, enstrophy 
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1 Introduction 



Determining the detailed properties of wall-bounded turbulent shear flows 
has been the focus of considerable experimental and computational research 
(see, e.g., [111213] for reviews). Although prior studies have provided insights 
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into low-order statistics and coherent structures in these flows, higher-order 
velocity gradient statistics remain relatively unexplored. Velocity gradients re- 
flect the structure and properties of the turbulent small scales, [I] and their 
higher-order statistics are particularly sensitive to the large-amplitude, inter- 
mittent fluctuations characteristic of high- Reynolds number flows. While such 
higher-order statistics have been studied in homogeneous isotropic turbulence 
using direct numerical simulations (DNS) for several decades, [5|6|7|8|9] it has 
only recently become feasible to carry out similar analyses in wall-bounded 
flows. This is due, in large part, to the substantial computational resources 
required to resolve the smallest scales of the turbulence. Progress in simulat- 
ing wall bounded flows has recently been made, however, in Ref. [10], where 
velocity gradient moments up to fourth order are examined using DNS of 
fully-developed turbulent channel flow at _Re T = 180, where Re T = u T L/u, u T is 
the friction velocity, L is the half-width of the channel, and v is the kinematic 
viscosity. 

The present paper refines and significantly extends the prior study in Ref. 
[TO] by using highly resolved turbulent channel flow DNS at Re T = 180 and 
Re T = 381 to examine both the height and Reynolds number dependence of 
higher-order velocity gradient statistics. Particular emphasis is placed on the 
moments of the energy dissipation rate 

e = 21/5^.55,. , (1) 



and the local enstrophy (or enstrophy density) 

n = \u>u> . (2) 

Taken together, these quantities characterize the straining and rotation asso- 
ciated with small-scale turbulent fluctuations. [8] The fluctuating strain rate 
tensor, S'^, in Eq. (CQ) is given in terms of the fluctuating velocity, v! si = Ui — (v,i) 
(where (■) is a z-dependent average over time and x-y planes parallel to the 
channel walls), as 



2 \dxj dxi 



S',, = -J^ + ^), (3) 



and the fluctuating vorticity, u'^ in Eq. (J2]) is given by 



where is the cyclic permutation tensor. The moments of e and Q, denoted 



2 



(e n ) and (fi n ), reflect properties of small-scale, high-amplitude fluctuations in 
the velocity gradient. For example, (e 4 ) / (e 2 ) 2 and (Q 4 ) / (Q 2 ) 2 can be used to 
quantitatively assess the degree of small-scale intermittency in the flow. [1] 

In order to measure moments of e and Q up to fourth order, the resolutions 
used in the present simulations are substantially higher than in prior sim- 
ulations at similar values of Re T . Recent studies of homogeneous isotropic 
turbulence [8lfT5] and turbulent channel flow [10] have shown that very fine 
resolutions are required to fully resolve higher-order statistics of e and Q. In 
particular, using simulations of a Re T = 180 channel flow at three different 
resolutions, it was shown in Ref. [10] that the lowest resolution was sufficient 
to measure low-order statistics such as mean velocities and Reynolds stresses, 
u'jii'y Finer resolutions, however, were required to accurately measure higher- 
order velocity gradient statistics such as the mean value of e, particularly 
near the channel walls. The resolutions used in the present simulations are 
even higher than those in Ref. [10], and allow calculation of up to fourth order 
moments of e and Q. 

A wide range of wall distances are examined in the present study in order to 
determine the height dependence of the statistics. Turbulent channel flows are 
typically divided into four regions; the viscous sublayer for z + <5, the buffer 
region for 5 < z + < 30, the logarithmic layer for z + > 30 and z < 0.3L, and 
the outer flow for z >0.3L, [TTJ where z + = zu T /v and z is the coordinate in 
the wall-normal direction. For both values of Re T considered here, statistics 
are examined from z + = 2, well within the viscous sublayer, to the channel 
centerline. The inhomogeneity in the z-direction requires that the statistics of 
e and Q be calculated in planes parallel to the channel walls, in order to allow 
averaging over homogeneous flow directions. Compared to studies of three- 
dimensional homogeneous isotropic turbulence, where full volume averaging 
is possible, this places significant restrictions on the statistical convergence of 
the results. Consequently, particular attention is paid in the following to the 
convergence of the moments of e and Q as a function of z + and Re T . 

While variations in the moments of e and Q with Re T and z + are important 
for understanding the small-scale structure of the channel flow, the presence 
of coherent vortical structures is also expected to play a role in determin- 
ing these moments. Theodorsen [17] first proposed a hairpin shape for these 
structures, and subsequent experimental (e.g., Refs. [T8|19|20] ) and numerical 
studies (e.g., Refs. [T2|21|22j ) have characterized hairpin vortices throughout 
wall bounded flows (see Ref. [2] for a review). Although these vortices have 
been identified as important in fluid transport and the generation of Reynolds 
stresses, [22] relatively little is known about their contribution to the higher- 
order, small-scale statistics of e and Q. In the following, we examine this issue 
using conditional analyses of the moments of e and Vt based on locations of 
intense rotation in the flow. We also consider the moments at locations away 
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from regions of intense rotation in order to determine the extent to which 
differences in the moments as a function of Re T and z + can be attributed to 
vortical structures. 

The manuscript is organized as follows. Details of the numerical simulations 
are presented in the next section. The moments of e and Q are then presented 
up to fourth order for both values of Re T , accompanied by an analysis of the 
statistical convergence of these moments. Conditional analyses based on loca- 
tions of intense rotation are outlined in Section HI The method by which these 
locations are identified is briefly discussed, and results from conditional analy- 
ses of the e and Q moments are presented. Finally, a summary and conclusions 
are provided at the end. 



2 Description of Numerical Simulations 



The numerical simulations used in the present study solve the incompressible 
Navier Stokes equations for a fully developed turbulent channel flow. These 
equations are written in non-dimensional form as 



^ = °' (5) 
duj | u dui _ dp ^ 1 d 2 Uj 
dt 3 dxj dxi Redxjdxj ' 

where itj is the velocity field and p is the kinematic pressure. The Reynolds 
number, Re, is given as Re = UL/u, where U is the total mean velocity in the 
channel. 

As described in Ref. [10], the Navier-Stokes equations in Eqs. (jSJ) and (jSJ) 
are decomposed in poloidal-toroidal form and then solved using a pseudo- 
spectral method. This method uses Fourier expansions in the x and y di- 
rections, which are parallel to the channel walls, and Chebyshev polynomial 
expansions in the z direction, normal to the walls. [231124] De-aliasing with the 
2/3 rule is applied in all three directions. The time discretization combines 
an implicit backward-differentiation method for linear terms with an explicit 
Adams-Bashforth method for nonlinear terms, and is second-order accurate. 
The method can be symbolically written as 

q fn+1 a fn I m— 1 

- — 2 [ t 1 = cr +l + 2Af(n - mr^ 1 ) , (7) 

where £ is a linear operator, M represents nonlinear terms, and At is the time 
step. The approximation for the time derivative, df /dt, on the left side of Eq. 
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Re T Re N x xN y xN z N t Ax + Ay + Az+ 
180 2800 512x512x1025 403 4.4 2.2 0.55 
381 6667 1024x512x1025 134 2.3 2.3 1.2 

Table 1 

Friction Reynolds number Re T , global Reynolds number Re = UL/is, grid dimen- 
sions N x x N y x N z , number of temporal snapshots, Nt, horizontal resolutions, 
Ax + = Axu T /u and Ay + = Ayu T /v, and vertical resolution at the centerline, 
Az£ = Az(z/L = l)u T /u, for the two simulations analyzed herein. 

(J7J) gives / at time level n + 1 using values at the two previous levels, denoted 
f n and f n ~ l . The simulations are parallelized using MPI. 

The channel flow simulation domain consists of a rectangular box of non- 
dimensional size (L x X Ly X L z ) / L = 47r x 27r x 2 for Re T = 180 and 2ir x ir x 2 
for Re T = 381. In Ref. [10], the Re T = 180 case was examined at a maximum 
resolution of N x xN y xN z = 512x512x513 collocation points in physical space. 
Due to de-aliasing, the number of Fourier or Chebyshev modes is 2/3 the 
grid numbers N XjVtZ . In the present study, this resolution has been increased 
to 512 x 512 x 1025 collocation points in order to achieve the best possible 
convergence of the velocity gradient statistics. The resolution in the _Re T = 381 
case, which was not examined in Ref. [10], is 1024x512x1025 collocation points 
(see Table [p. These parameters give grid spacings that are approximately 
two to five times finer than in prior simulations at comparable values of Re T . 
The Re T = 180 channel flow simulated by Kim et al. [12], for example, used 
192 x 160 x 129 grid points in a domain of size 47rx27rx2 (where, consistent 
with the present study, z is used as the wall-normal coordinate). Similarly, the 
Re T = 395 channel simulation by Moser et al. [13]) was performed on a grid 
consisting of 256 x 192 x 193 points in a domain of size 2ir x 7r x 2. 

Temporal snapshots of the flow field are stored every 0.1 convective time units, 
L/U, for the Re r = 180 simulation, with a total of N t = 403 snapshots. The 
time interval between snapshots for the _Re T = 381 simulation is 0.19 convec- 
tive time units, with Nt = 134 total snapshots. The statistical analysis for a 
single horizontal plane (and its symmetry counterpart in the upper half of the 
channel) can thus be made over a set of up to 2.1 x 10 s data points. Due to 
the larger N t for the Re T = 180 case and the greater N x for Re T = 381 (see 
Tabled]), we use a similar amount of data in the analysis of each value of Re T . 

Figure Q] shows the mean velocity in the x direction and the Reynolds shear 
stress Ti 3 = u' x u f z + = u' x u' z /ul as a function of z + for i?e r = 180 at a resolution 
of 128x128x129. Even for this relatively low resolution, the simulation results 
are in good agreement with recent channel flow simulations by del Alamo and 
Jimenez [16]. The high resolution of the present simulations further allows the 
higher-order statistics of e and Q to be calculated. 
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Fig. 1. Mean streamwise velocity u + (a) and Reynolds shear stress t^ = u' x u' z (b) 
from i?e T = 180 numerical simulation at resolution 128x128x129. Results show good 
agreement with DNS of fully-developed turbulent channel flow by del Alamo and 
Jimenez [TB] for i?e r = 180. 

3 Statistics of Energy Dissipation Rate and Local Enstrophy 

The statistics of e and Q are examined in the following at various wall dis- 
tances, z + , by carrying out the analysis in x-y planes parallel to the channel 
walls. Within each of these planes, the flow is essentially homogeneous, thus 
allowing an examination of the flow statistics similar to that employed for 
homogeneous isotropic turbulence. In the following, the average (■) denotes 
an x-y average at a particular value of z + . It is always combined with an 
arithmetic average over the full sequence of temporal snapshots. Statistical 
convergence is further improved by using symmetric planes from both the top 
and bottom halves of the channel, where the velocity and velocity gradient 
fields in the top half are reflected about the centerline. 

3. 1 Distributions of the Dissipation and Enstropy 

Figure [2] shows probability density functions (pdfs) of e/(e) and for 
both values of Re T at wall distances spanning the viscous sublayer (z + = 2) to 
the channel centerline, where the averages (e) and (Q) are functions of z + . The 
pdfs of both e and Q vary as the wall is approached from the centerline, and 
Figure [2] shows that the degree of intermittency, as indicated by the width of 
the pdf tails and the corresponding probability of obtaining large amplitudes 
of e and Q, generally decreases as the wall is approached, with a minimum in 
the buffer layer at approximately z + = 10. Immediately at the wall, for z + = 2, 
however, the pdfs of e and Q for both values of Re T are wider than those for 
z + = 10. For Re T = 381, the tails of the e and Q pdfs are most prominent in 
the logarithmic layer for z + ?y 60 — 90, while for Re T = 180 the most prominent 
tails occur further from the wall at z + ~90 — 150. This is, in turn, indicative 
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Fig. 2. Probability distribution functions of dissipation rate s/{e) and local enstro- 
phy Q/(n) for i?e r = 180 (e: (a) and ft: (b)) and #e T = 381 (e: (c) and ft: (d)). 

of increased intermittency at these values of z + , and may be connected to the 
bursting of coherent vortical structures away from the wall. [2] 

With respect to the differences between the pdfs of e and Q, Figure [2] shows 
that for both values of Re T and wall distances down to z + ~40, the tails of the 
Q pdfs are more pronounced than those for the e pdfs. This is consistent with 
prior results in homogeneous isotropic turbulence j8] and turbulent channel 
flow [ID] . For z + <40, however, the pdfs of e and Q are more similar, indicating 
a closer correspondence between the statistics of e and Q near the wall. 

The connection between e and Q at each value of z + can be further examined 
using joint probability distributions of e and Q, shown in Figure [3] for Re T = 
381. Similar joint pdfs are shown for Re T = 180 in Ref. [10]. The joint pdfs, 
denoted P(e,Q), are normalized by P(e) and P(£l) as [TO] 



n(e,fi) 



P(e)P(n) ' 



where values of II (e, f2) greater than 1 indicate a higher correlation between e 
and f2 than if the two quantities were statistically independent. [25] Consistent 
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Fig. 3. Joint pdfs of s and f2 at six values of z + for Re T = 381. The joint pdfs are 
normalized by the one-dimensional pdfs of e and Q, as in Eq. ([8]), and log 10 [n(e, f2)] 
is shown. The color contours are the same in all panels, and range from 10 -4 (blue) 
to 10 3 (red). The black diagonal line in panel (f) corresponds to the relation e=2v£l 
obtained from Eq. ([9]). 



with prior results for Re T = 180 [10] and with simulations of homogeneous 
isotropic turbulence, [25] Figure |3] shows that intense e and Q are correlated 
at all values of z + . As the wall is approached, however, Figure [3] shows that 
the support of the joint pdfs decreases. As noted in Ref. [10] . the averages of 
e and fl are connected in a wall-bounded shear flow by 



d 2 (u' 2 ) 

(e) = 2u{Q) + 2v-^l 



(9) 



Since the second term on the right-side of (Q is small compared to 2u(Q) in 
the channel flow, [TUJ we obtain (e) «2i/(fi). Figure [3] shows that the values of 
e and Vt fall increasingly close to this relation as the wall is approached. These 
results at Re T = 381 are qualitatively consistent with the findings in Ref. [10] 
for Re T = 180. 



3.2 Moments of the Dissipation and Enstrophy 



The statistical convergence of a moment ((p n ) can be assessed from plots of 
(p n P(ip) versus </?, where </? is an arbitrary velocity gradient quantity such as e 
or Q. [8] These distributions are shown for n = 3 and n = 4 for the _Re r = 180 
and Re T = 381 simulations in Figures @] and El respectively. Figure H] shows 
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Fig. 4. Moment convergence results for e and SI for i?e r = 180. The n = 3 and n = 4 
moments are shown for e (a) and (c), and (b) and (d), respectively. 

that, for Re T = 180, the n = 3 moments are converged for both e and at 
all z + , as indicated by the decrease in the curves to zero for large e/(e) and 
0/ (f2). While the n = 4 moments of e are also converged at all z + , the n = 4 
moments of Q show a lack of convergence for z + « 90—150. These wall distances 
correspond to the pdfs with the longest tails in Figure [2(b), and the lack of 
convergence shown in Figure 11(d) is likely due to the bursting of vortical 
structures across the channel. This would also explain why the higher-order 
Q moments are more strongly affected than the e moments, and suggests that 
the non-converged fourth moments of Q in Figure 11(d) may have a physical 
origin. 

For i?e T = 381, Figures [5(a) and (c) show that the moments of e are adequately 
resolved up to n = 4 at all values of z + . Figures [5(b) and (d) further show that 
the n = 3 and n = 4 moments of Q are generally better converged than those 
for i?e r = 180 in Figure HJ Given the fact that approximately the same amount 
of data is used in the analysis of both runs, this suggests that the cross- 
channel bursting of vortex structures may be less prevalent for higher Re T . In 
particular, there is an extended bulk region for the higher Reynolds number 
case that separates the two logarithmic layers in the upper and lower halves 
of the channel. 
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Fig. 5. Moment convergence results for e and SI for Re T = 381. The n = 3 and n = 4 
moments are shown for e (a) and (c), and (b) and (d), respectively. 

Figure |6] shows the resulting moments for n = 2— 4 for the Re T = 180 and Re T = 
381 simulations throughout the entire channel. Estimates of 90% confidence 
intervals for the calculated moments have been included in Figure |6] using a 
moving-block bootstrap method [26] on time-series of (s n ) xy}S and {fl n ) xytS , 
where (-)xy,s denotes a time- and ^-dependent average over x-y planes and 
symmetric halves of the channel (s). The block sizes used in the bootstrap 
analyses are approximately equal to one convective time unit (L/U) for both 
values of Re T , giving ten snapshots per block for Re T = 180 and five snapshots 
per block for i?e T = 381. 

For both values of Re T , Figure |6] shows that the e and Q moments are large very 
close to the wall at z + = 2, but reach their minimum values at z + = 10 within 
the buffer layer. The large amplitudes of the moments at z + = 2 arise from a 
combination of effects, including normalization of the moments by the mean 
values of e and Q, and the fact that z + = 2 corresponds to locations within the 
viscous sublayer where the flow is not fully turbulent. Moreover, while velocity 
fluctuations become small near the wall, fluctuating velocity gradients become 
very steep, resulting in very large values of e and Q at z + = 2. For both values 
of Re T , the moments of e and Q are similar up to z + ^30 — 40, but for larger 
values of z + , the moments of Q are substantially greater than those for e. This 
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Fig. 6. Moments (e n )/(e) n (black lines) and (Q n )/(Q,) n (red lines) as a function 
of z + for n = 2 — 4 in the Re T = 180 (a)-(c) and Re T = 381 (d)-(f) simulations. 
Legend is shown in panel (c). Error bars show estimates of 90% confidence intervals 
obtained using a moving-block bootstrap method, [26] with block sizes equal to one 
convective time unit (L/U) for both values of Re T . 

last result is consistent with the wider pdfs of Q compared to those for e for 
z + > 40 in Figure [2j and has also been observed in studies of homogeneous 
isotropic turbulence. [8] 

While Figure |6] shows that the dependence of the moments on z + is similar 
in certain respects for Re T = 180 and Re T = 381, there are slight differences 
between the two values of Re T . For i?e r = 381, the moments of e reach a local 
maximum near z + m 50 before decreasing until z + ~ 100. For larger values of 
2 + , the e moments remain relatively constant to the channel flow centerline. 
For Re T = 180, by contrast, the moments of e do not show a pronounced 
peak anywhere in the channel and remain relatively constant outside of the 
buffer layer. The local maxima in the e moments at z + ~ 50 for Re T = 381 
could indicate a qualitative change in the small-scale turbulence, in accordance 
with Yakhot et al, [27] where it was suggested that small-scale fluctuations 
outside the buffer layer should be similar to those in isotropic turbulence. The 
invariance of the e and Q moments for z + > 100 in the Re T = 381 simulation 
provides support for this idea. 

Figures E)^a)-(c) further show that the moments of Q for Re T = 180 reach a 
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maximum near z + m 110— 150 before decreasing again at the channel centerline. 
This pronounced maximum is likely due to the bursting of coherent vortices, 
but is not observed for Re T = 381 in Figures E^d)-(f), where the moments of Q 
remain relatively constant for z + >100. The error bars in Figure 0(a)- (c) for 
Re T = 180 are, however, largest in this range of z + , indicating that, particularly 
for the fourth moments, additional data is required to fully capture the strong 
temporal variability of the vorticity field created by these bursting structures. 



4 Conditional Analyses of Dissipation and Enstrophy Moments 

The statistics of e and Q in the previous section are obtained using all points 
in each plane of the channel. Certain features of these statistics, such as the 
local maxima in the moments of Q for Re T = 180 at z + ~ 110 — 150 and the 
differences between the moments for e and Q outside of the buffer layer, may 
be due, in part, to the presence of intense vortical structures in the channel. 
This issue can be examined by calculating the statistical moments at points 
both within and outside regions of strong rotation in the flow, which are taken 
here to correspond to vortical structure locations. 



4-1 Identification of Intense Vortical Structures 

There has been considerable research over the last several decades on the most 
appropriate methods by which to identify intense vortical structures in turbu- 
lent flows. Chakraborty et al. [28] showed, however, that most vortical struc- 
ture identification procedures give similar results for homogeneous isotropic 
turbulence. Since our primary interest here is not in characterizing the prop- 
erties of the structures or assessing the merits of various identification pro- 
cedures, we simply classify vortical structures as intense, rotation-dominated 
regions of the flow. The analysis of the e and Q moments is then carried out 
only over points that fall within these regions, or points completely outside 
these regions, which we term the background flow. 

As first proposed by Hunt et al, [29] rotation-dominated regions can be iden- 
tified as locations where Q>0, where Q is the second invariant of the velocity 
gradient tensor, Aij = dui/dxj, and is written for an incompressible flow as 




(10) 
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The strain rate tensor, Sij, is given by 



and is the anti-symmetric part of Ay given by 

(12) 



Here we use a closely related identification method, originally proposed by 
Zhou ef al, (30] that requires Ajj to have two complex conjugate eigenvalues, 
\ r ±i\ ci , and that the complex part, A C j>0, be larger than a prescribed cutoff. 
The restriction that Aij have two complex conjugate eigenvalues is equivalent 
to requiring that A>0, where A is given by [22] 

A = Q 3 + ^R 2 , (13) 



and 



R — —— (SijSjkSki + 3RijRjkSki) , (14) 



is the third invariant of Ay for an incompressible flow. Regions of intense 
rotation - and, by extension, vortical structure locations - are then identified 
by requiring that \ C i>a (Ki) max , where (Ki) max is the maximum value of A C j 
at each value of z + . The pre-factor a determines, in large part, the number 
and magnitude of vortical structures identified in the flow, and we use a = 
0.05 herein. The resulting locations of intense rotation, identified as points 
where A C j is large, are generally similar to those obtained using other criteria, 
including the Q criterion. 



Due to the presence of a mean shear in the channel, the entire vortical structure 
identification procedure is carried out using the fluctuating velocity gradient, 
A\j = du\l 'dxj. This approach has been used previously by Robinson [31] and 
Pirozzoli et al, [22] and is particularly important near the channel walls where 
the mean shear can generate large vorticity not associated with coherent vor- 
tical structures. The fluctuating tensors S'^, given in Eq. (j3J), and R'^, given 
by 



(15) 
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Fig. 7. Schematic showing vorticity orientation angles defined in Eq. (I16p , 

are thus used in the expressions for Q, R, and A in Eqs. (TTU]) . ( IT41 . and (fTBl . 
The difference between using the full velocity gradient tensor and the fluctu- 
ating tensor is small for much of the channel, and only becomes significant in 
the near- wall region where the mean shear is large. 

4-2 Orientation of Intense Vortical Structures 

In the present analysis, the orientation of intense vortical structures is inferred 
from the orientation of the fluctuating vorticity field, cu-, at rotation-dominated 
locations in the flow. This is in contrast to prior approaches (e.g., Ref. [20] ) 
which have attempted to treat coherent vortices as connected regions with a 
single orientation associated with the structure as a whole. Following Pirozzoli 
et al, [22] the orientation of the vorticity is characterized using the angles 

% = tan" 1 (^j , 9 e = sin" 1 ^ , (16) 

where tu' = (u;-^) 1 ^ 2 , which are shown schematically in Figure [71 The orien- 
tation of intense vortical structures can then be determined from the joint 
pdfs of 9 xy and 9 e , denoted P(9 xy ,9 e ). Since an isotropic vorticity field gives 
P(9 xy ,9 e ) ~cos(6' e ), we consider normalized joint pdfs of P(9 xy ,9 e )/ cos(9 e ), 
following the approach used in Ref. [22j . 

Joint pdfs of P(9 xy , 9 e )/ cos(6 l e ) for Re T = 381 in Figure [H] show that the dis- 
tribution of vorticity within rotation-dominated regions varies substantially 
as the wall is approached from the channel centerline. There is essentially no 
preferred orientation of the vorticity near the centerline, as shown by the lack 
of any clear maxima or minima in Figure [8]^a). At z + m 120 in Figure EJ^b), 
however, the vorticity shows a weak preference for # e ^±45° and #^^±90°, 
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Fig. 8. Joint pdfs of the elevation angle 9 e and the orientation angle 6 xy , defined in 



Eq. (I16p . at six values of z + for Re T = 381. Only vorticity identified to be within 
vortical structures is included in the pdfs (i.e. points for which A c j > 0-05(X c i) max , 
where (X c i)max is determined at each z + ). Contours of log 10 [P{0 xy ,0 e )/ cos# e ] are 
plotted, with levels from 10 -3 (blue) to 10° (red). 

approximately corresponding to the "necks" of hairpin vortices discussed in 
previous studies. [20)122] As the wall is approached, Figure shows that 
the peaks in the distributions begin to shift towards 9 e = 0°, corresponding to 
quasi-streamwise vortices [TTf2"2"] which are oriented in the x-direction parallel 
to the channel walls. At the same time, the probability for 9 xy > 1 90° | also be- 
gins to increase, and by z + = 10 in Figure [8]^e) there is a strong probability of 
obtaining vorticity with 9 xy close to ±180°. At z + = 2 in Figure [HJ^f), much of 
the vorticity is oriented in the spanwise direction parallel to the channel walls 



(giving 9 e 



9 « 

U X y ~ 



0°). There are also less pronounced peaks at \9 xy \ = 180°, 



which are due to weaker vorticity not associated with intense structures (as 
indicated by the disappearance of these secondary peaks when using larger 
values of a). There is thus a rapid change in the orientation of intense vortices 
between the buffer layer and viscous sublayer, a result that is mirrored in the 
large changes in the moments of e and Q between z + = 10 and z + = 2 in Figure 

El 



4-3 Conditioned Moments of the Dissipation and Enstrophy 



Figure M shows the moments of e and Q conditioned on vortical structure 
locations, which are identified as rotation-dominated points in the flow using 
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Fig. 9. Moments (e n )/{e) n and (n n )/(Q) n for n = 2-4 in the Re T = 180 (a)-(c) and 
Re T = 381 (d)-(f) simulations. The curves correspond to moments obtained from 
the full fields (black lines, denoted "Full"), points at rotation-dominated locations 
(red lines, denoted "Vortex"), and points where the rotation is weak (blue lines, 
denoted "NV"). 

the procedure outlined in Section 14.11 Moments are also shown for points in 
the background field, where the rotation is small (see the A C j criterion defined 
immediately below Eq. ( jl~4l ). The moments for both subsets are normalized 
with respect to the means of e and Q obtained from the full field. 

For both e and Q, Figure M shows that the moments obtained at locations of 
intense rotation exceed those for both the full and background fields. Although 
the increase in the Q moments at these locations is to be expected (since the 
identification procedure essentially selects points with large vorticity magni- 
tude), the accompanying increase in e can be understood from the joint pdfs 
in Figure [31 These pdfs show that intense e and Q are statistically correlated, 
and thus the moments of e at locations of intense rotation tend to be larger 
than the corresponding full field values. From studies of isotropic turbulence 
(see, e.g., [32]), it is also well-known that large amplitudes of both e and fl 
appear in close spatial proximity. At the same time, however, it should be 
noted that the moments of e remain smaller than those of Q at all z + , for all 
orders, and at both values of Re T . 

Substantial differences between the moments of Q and e are observed in the 
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background field, where the rotation is small. In particular, there is near agree- 
ment between the e moments in the full and background fields. This indicates 
that locations of intense rotation, and hence vortical structures, do not con- 
tribute substantially to the full statistics of e. By contrast, Figure M shows 
that the background moments of fl are significantly smaller than the corre- 
sponding full field values, indicating that locations within vortical structures 
make a substantial contribution to the full field moments of fl. Contrary to 
the moments in the full field and at rotation dominated locations, the local 
maxima in the moments of fl at z + ~ 150 for Re T = 180 are absent in the 
background field. This reinforces the connection between these maxima and 
intense vortices in the flow. Such pronounced maxima are not observed in any 
of the fields shown in Figure [9] for Re T = 381, and it remains to be seen in 
future studies whether these results undergo additional changes for even larger 
Re T . 



5 Summary and Conclusions 

The present high-resolution DNS study of fully-developed turbulent channel 
flow has examined velocity gradient statistics as a function of wall distance, 
z + , at friction Reynolds numbers Re T = 180 and _Re r = 381. An emphasis has 
been placed on the statistics, and in particular the higher-order moments, of 
the energy dissipation rate, e, and the local enstrophy, fl. The probability 
density functions (pdfs) of e and fl for the Re T = 381 case are qualitatively 
similar to previous results obtained for Re T = 180. [10] The pdfs of both e 
and fl generally become less intermittent as the wall is approached from the 
channel centerline, as indicated by the less broad tails of the pdfs. For z + >40, 
the pdfs of fl have wider tails than the pdfs of e, although the pdfs for both 
quantities are similar near the channel walls. Joint pdfs show that both large 
and small values of e and fl are correlated, and that the support of the pdfs 
is reduced as the wall is approached. 

The high resolution of the present simulations has allowed moments of e and 
fl up to fourth order to be calculated. The moments of e and fl are similar in 
the viscous sublayer and buffer layer, but, for locations further from the wall, 
the moments of fl are substantially larger than those of e. Reynolds number 
effects are observed in the moments of both e and fl. In particular, there is 
a local maximum in the moments of e at the beginning of the logarithmic 
layer for Re T = 381, which is not present in the moments of e for Re T = 180. 
For both values of Re T , the moments of e remain relatively constant from 
z + 100 to the channel centerline. The fl moments also remain relatively 
constant over this range for Re T = 381, but for Re T = 180 the moments of fl 
increase and reach a maximum at z + ps 110 — 150, before decreasing at the 
channel centerline. The maxima in the moments of fl for the smallest Re T are 
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most likely due to bursting of vortical structures across the channel, which is 
an effect due to the low Reynolds number. In particular, it is possible that both 
wall regions may not be fully decoupled. This bursting causes large temporal 
variations of the local enstrophy, which complicates the statistical convergence 
of higher order moments. This lack of convergence may thus have a physical 
fingerprint, and substantially more temporal snapshots (requiring significant 
additional computational effort) are necessary to resolve this issue in future 
investigations. 

Using conditional analyses based on regions of intense rotation, which are 
taken here to correspond to vortical structure locations, the moments of both 
e and Q are shown to be larger in rotation-dominated regions than in the 
full field. The increase in e, in particular, is due to the correlation between 
events of intense e and Q shown in Figure [3J At the same time, moments 
calculated in the background field, where the rotation is not intense, indicate 
that vortical structures make only a small contribution to the moments of e in 
the full field. Differences between the moments for the two values of Re T in the 
conditional intense rotation fields are similar to the differences in the full field. 
In particular, the conditional moments of e and Q are relatively constant for 
z + > 100 in the Re T = 381 case, but local maxima, particularly in the moments 
of Q, are still observed at z + f» 150 for the Re T = 180 case. These maxima are 
not present, however, in the Q moments in the background field for Re T = 180. 

The next step in the analysis is to extend the study of the Reynolds num- 
ber dependence of the statistics, which requires further high-resolution data 
at larger values of Re T . This will allow us to make comparisons with similar 
studies in isotropic turbulence, (15] and to identify quantitative differences in 
the statistics close to the channel walls. These simulations, which are currently 
in progress, require substantial computational resources and will be reported 
in the near future. 
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